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Abstract. Taking up a method devised by Taylor and Wheeler and collecting 
pieces of their work we offer a self-contained derivation of the formulae giving both 
the precession of the orbit of a planet around the Sun and the deflection angle of 
a light pulse passing near the Sun in the framework of General Relativity. The 
demonstration uses only elementary algebra without resorting to tensor formalism. 
No prior knowledge in relativity is needed to follow the presentation. 
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1. Introduction 

The deviation of light rays near the Sun and the cliange in the orientation of Mercury's 
orbit with each revolution were the first dramatic confirmations of the correctness of 
Einstein's theory of general relativity. Surprisingly it is uneasy to find the detailed 
calculations of those phenomena in the specialized scientific literature. For example 
even the comprehensive reference book Gravitation by Misner et al (1973) does not 
provide the reader with full explicit demonstrations of the effects. Truly it would seem 
that there is no simple way to conduct the complete calculation since it is necessary to 
dive into the full theory of general relativity. Fortunately, in their book about black 
holes and General Relativity, Taylor and Wheeler (2000) have succeeded in deriving 
the needed equations of motion without appealing to tensor formalism. To determine 
the motion of a free object in curved geometry, they use what they call a "Principle 
of Extremal Aging" along the path of the particle. That principle directly leads to the 
desired equations while using only simple algebra. 

Exploiting their astute technique and collecting together parts of their work we 
carry out a complete presentation of the legendary effects predicted by Einstein when 
founding General Relativity, namely the deviation of stellar light near the sun's mass 
and the advance of the perihelion of Mercury. The aim of the present paper is to supply 
teachers, students and beginners in General Relativity with a document both yielding in 
a single place a full account of those relativistic phenomena and still remaining simple 
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and accessible in spite of the difficulty of the problem. In fact in a sole article we begin 
with the derivation of the equations of motion in the curved spacetime surrounding 
a massive object and we end with the numerical evaluation of the magnitude of the 
physical effects under consideration. The calculation is thus entirely self-contained. 
Now, as far as we are aware, such a kind of document does not exist yet and thus deserves 
being published, even if it does not offer new material, as almost all the demonstrations 
shown here may be found in Taylor and Wheeler's book. 

We first closely follow Taylor and Wheeler's treatment in estabhshing the required 
formulae of motion in Schwarzschild geometry. But we depart from their approach when 
calculating the advance of Mercury's perihelion. Carrying out a calculation adapted 
from Wald (1984), Taylor and Wheeler analyze the orbit of the planet by splitting 
its periodic motion into two components: an angular orbit and a radial oscillation. 
While the angular motion is described as nearly circular, the small inward and outward 
displacement is considered as a simple harmonic oscillator moving radially in a potential 
well. In this context it is found that the two corresponding periods differ by a tiny 
amount. As the oscillation takes place more slowly, it becomes late with respect to the 
period of revolution and this gives rise to the precession phenomenon. Contrary to this 
strategy we prefer to proceed more directly by adopting the algebraic method used by 
Weinberg (1972) in his book on General Relativity. The corresponding procedure allows 
the orbit to be elliptical and docs not require the development of additional concepts. 
Actually, as they arc written in polar coordinates, the equations of motion immediately 
yield the progression of the azimuthal angle as a function of the radial coordinate r 
and thus make it possible to evaluate the change in this angle over one revolution. In 
particular the formulae show that from one passage at the perihelion to the next one the 
azimuthal angle increases by a little more than 27r, which means that the perihelion 
does advance along the orbit. 

Additionally, taking advantage of the expressions at hand we rediscover Kepler's 
third law relating the period of revolution of a planet to the length of its semi-major 
axis. 

We finally write the formulae determining the trajectory of a photon around a 
massive body and compute the change in the direction of a starlight pulse passing near 
the Sun. The integration of the equations readily shows that the azimuthal angle of the 
photon increases by more than tt along its path (tt corresponding to a straight line in 
the absence of mass), the additional quantity representing the relativistic deflection. 

Incidentally, as the paper is intended to be fully accessible to students and (again) 
self-contained, some elementary notions of physics, especially in relativity, are recalled 
here and there. As a result no prerequisite in relativity is needed to study this article. 

2. Principle of the calculation 

The orbit of Mercury, the closest planet to the Sun, was known as early as 1859 to 
change with time in a partly unexplained way. It was observed that the perihelion (the 
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orbital point closest to the Sun) moves a little around the Sun, advancing at a tiny rate 
with each orbit. But whereas Newtonian mechanics was able to account for the major 
part of this change, a residual of 43 seconds of arc per century remained unexplained. 
It is this residual that was correctly predicted by Einstein's theory. Likewise Arthur 
Eddington's verification of the deflection of starlight near the Sun in the total solar 
eclipse of 29th May, 1919, gave Albert Einstein an immediate notoriety. Now while 
the quahtative aspect and the historical impact of both the precession of Mercury's 
perihelion and the deviation of light rays near the Sun are widely (and legitimately) 
described in the scientific literature as proofs of the rightness of General Relativity, the 
detailed calculations of the magnitude of those tiny relativistic effects are rarely given. 
The present paper is consecrated to this "numerical exercise" . 

General relativity teaches us that mass particles or photons propagating freely 
through space follow what is called a "geodesic" of spacetime. Thus we have to examine 
the following points: 

• what is a "geodesic"? 

• which equations govern a geodesic? 

• find the solutions of those equations for a mass particle and determine the orbit of 
a planet (which is actually an ellipse in a first approximation: thank you Newton!); 

• compute the tiny angle of rotation of the ellipse after the completion of one orbit; 

• by letting the mass of the particle tend towards zero, find the trajectory of light 
near the Sun; 

• compute the angle between the initial and final directions of a light ray grazing the 
solar disk. 

3. Metrics for fiat spacetime in special relativity 

Newtonian mechanics describes the motion of a particle in an absolute space with respect 
to an absolute time. The position of the moving object is located by its coordinates 
with respect to some frame of reference and is given as a function of time t. Special 
relativity affirms that there exists no absolute time and that time cannot be dissociated 
from space. The theory bases its reasoning on events, each event being characterized by 
a point M (where it happens!) and a time t (when it happens!). The events attached 
to a moving particle constitute what is called a worldline. 

Let us consider for example a spaceship moving freely through space, which means 
that all its motors are turned off. Let us imagine that regular flashes are emitted in 
accordance with a clock located inside the rocket and ticking at a constant rhythm. The 
time interval between two successive flashes will be denoted by r (this quantity is thus 
measured with respect to the proper time of the spaceship). Think now of another frame 
of reference as constituted by an ensemble of space beacons, also free from acceleration, 
at constant mutual distances from one another (each free beacon stays at the same 
distance from its neighbours). Every signal bears the indication of its position in space 
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(for instance by showing its distance from some point of origin) and holds its own clock. 
The clocks of this second frame are synchronized between them. Then in that latter 
frame the interval between two flashes (i.e. two events) is characterized by two numbers: 
the space interval r and the time interval t. To determine those two quantities it suffices 
to record which beacon faces flash j^l and which beacon faces flash while noting the 
times of those events as shown by the clocks attached to the beacons. 

Special relativity is based on the following principle. The proper time interval r 
between event #1 and event #2 is given by the formula 



and this quantity does not depend on the frame in which it is evaluated. In other words 
all observers agree on the value of r computed by formula (1), although the values of r 
and t differ from one system of reference to another. 

Unless otherwise indicated, distances will be measured in units of time, as is 
often done in astronomy. We have chosen to do so in writing (1). If however 
distances r are expressed in conventional units, for instance in centimetres, then one 
should pass from the latter to our distance r expressed in seconds via the formula 
r(in seconds) = r(in centimetres) /c where c is the speed of light in conventional units, 
namely 3 x 10^° cm/ s. (Expressing distances and times with the same unit would amount 
to taking the speed of light equal to unity.) 

4. Metrics for curved spacetime around a rucissive object 

In general relativity the property of invariance of the proper interval with respect to 
a change of the coordinates remains valid but only locally, i.e. under the condition 
of staying in a sufficiently small region of spacetime (its size depends on the accuracy 
of the measurements). The new feature concerns the expression of the proper time 
previously given by (1). The coefficients entering this formula depend now on the 
point of spacetime under consideration and the resulting expression takes the name of 
metrics. In fact, and this is wonderful, the whole structure of spacetime, and especially 
its curvature, is included in the local expression of r and in the form of its coefficients. 

We are interested here in the structure of spacetime around the Sun. In order to 
describe the physics locally we consider two nearby events separated by infinitesimal 
amounts of time and space coordinates dt, dx, dy and dz. If space were flat, the metrics 
would have the form 



By working in spherical coordinates, in a plane containing the center of the Sun (this 
choice removes one spatial coordinate), that formula becomes 



T' = t'-r' 



(1) 



{drf = {dtf - {dxf - {dyf - {dzf 
which is usually written (and a little sloppily) by convention as 
= dt^ - dx^ - dy^ - dz'^ . 



(2) 



dr^ = dt^ - dr^ - r^d(p'^ , 
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where r denotes the distance to the center and (f) an azimuthal angle in the plane of the 
orbit. 

But spacetime around a center of attraction of mass M (for instance a black hole 
or the vicinity of the Sun) is not fiat. It is characterized by the Schwarzschild metric 
(see for instance Misner et al, 1973) 

dr^ = (1 - 2M/r)dt^ - (1 - 2M/r)-^dr^ - r^dcj)'^ . (3) 

The story is really fantastic: the whole structure of spacetime is embodied in 
this "simple" formula (3). Even the famous black hole lurks behind those apparently 
innocuous symbols. 

One question: in which units is the mass M expressed in that formula? It is seen 
that M has the dimension of a length, a quantity that we measure here in seconds. 
Therefore M will also be measured in seconds. The formula allowing to transform 
grams in seconds is 

M(in seconds) = (G'/c^)M(in grams 

where (G/c^) = 2, 5 x 10"^^ s/g. Thus the mass of the Sun is 2 x 10^^ g or 5 x 10"'' s 
or 1.5 km. 

5. The equations of a geodesic 

The metric, that is the formula expressing at a given point of spacetime the temporal 
interval between two nearby events, reveals the presence of curvature as soon as the 
expression deviates from formula (2) that corresponds to fiat euclidian space. This 
metric will allow us to find the properties of the motion of a test particle free from 
acceleration. Actually both special and general relativity teach us that between two 
given events Ei and E2 a freely moving body follows the path for which the time interval 
r is the maximum. By defining a geodesic as the path maximising the time interval, it 
can also be stated that a freely moving particle follows a geodesic of spacetime. 

Definition of a geodesic: the geodesic between two events Ei and E2 is the 

wordline for which the interval of proper time between Ei and E2 is maximum. 

That property of maximising the proper time will allow us to derive the equations 
of a geodesic. It will also yield the expressions of the energy and angular momentum of 
a particle in orbit around the centre of attraction. 

6. Energy of the peirticle 

Let us apply the principle of maximisation of the proper time interval in the following 
manner. Suppose that a free spatial ship (whose rockets are turned off) falls radially, 
therefore along a straight path, towards the central attractive mass. Imagine that 
three successive fiashes, with nearby time and space coordinates, are emitted inside the 
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spaceship. We observe those three events in some external frame. In that latter frame 
the event Ei consists in the emission of a flash at time t = when the spatial engine 
is located at radius ri. The flash E2 is emitted at time t when the cabin is at radius 
r2. The flash is emitted at time T when the cabin is at radius r^. The quantity T 
is assumed to be small. We then imagine that we vary the intermediate coordinates of 
E2. The principle of maximal aging says that the geodesic starting from Ei and ending 
at £'3 will pass through event E2 such that the proper time interval 

t = ta + tb, (4) 

is maximum. Here ta measures the interval over the first spacetime segment A, which 
connects Ei to E2 and tb measures the time interval over the second segment B, which 
connects E2 to E3. 

In order to avoid varying all quantities at the same time, we assume in this 
experiment that the locations of the radii ri, r2 and are fixed and that only the 
time t, at which the second flash is emitted, is allowed to change. According to (3) the 
interval of proper time over the first segment A is given by its square 

TA^ = (1 - 2M/rA)t^ + (terms without t) , (5) 

from which we deduce 

TAdTA = (1 - 2M/rA)tdt . (6) 

The lapse of time over segment B between the events E2 and £'3 is {T — t), and 
therefore the proper time duration tb is given by 

tb^ = (1 - 2M/r)(T -tf + (terms without t) , (7) 

from which we deduce 

T^dTB = -(1 - 2M/rB){T - t)dt . (8) 

To make the total time interval t = ta + tb maximum with respect to a variation 
dt of the time t, we write 

^ = ^ + ^ ^ f9) 
dt dt dt ' ^ ' 

Deducing drA and dTB from (6) and (8) and letting quite naturally t — tA and T—t — tB, 

we easily get 

(1 - 2M/rA)itA/rA) = (1 - 2M/rB)(tB/rB) . (10) 

The left side of that equation depends only on parameters characterizing the first 
segment A (which connects Ei to E2). The right side depends only on parameters 
related to the second segment B (which connects E2 to £3). 

We have discovered in equation (10) a quantity that is the same for both segment. 
This quantity is thus a constant of the motion for the free particle under consideration. 
For good physical reasons (especially to recover the formulae of special relativity), one 
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is led to identify that constant of motion as the ratio of the energy of the particle to its 
mass. We write this very important result under the form 

E/m^(l-2M/r)(dt/dT), (11) 

an expression in which we have returned to the differential notation for the intervals t 
and T. 

Incidentally we may notice that with the units we have chosen, energy E and mass 
M are expressed in the same unit (for instance the second or the centimetre). 

7. Angular momentum of the particle 

We have applied the principle of maximising the proper time interval by varying the 
time of the intermediate event E2. We now perform the same operation but this time 
we vary the angle (f) of that intermediate event. We recall that (f) measures the direction 
of the moving particle with respect to a direction chosen as the origin. We call it the 
azimuth. 

We again consider the three events consisting in the emission of flashes inside a 
spaceship floating freely in space. The flrst segment A connects event Ei to event E2. 
The second segment B connects E2 to £"3. The azimuthal angle of the flrst event is 
fixed at = 0. The angle of the last one is fixed at = The intermediate azimuth 
is taken as the variable (p. Again in order not to vary everything at the same time, we 
assume that the radius r at which the second flash is emitted stays constant. 

We follow the same chain of reasoning as in the previous section. Prom the 
metric (3), the time interval ta over the first segment is given by its square 

ta^ = -ta^^)^ + (terms without 0) , (12) 
and the interval tb over the second by 

tb^ = -rij^($ - 0)^ + (terms without 0) , (13) 
from which we get 

TAdTA^ -rA^(t>d(t) (14) 
TBdTB^rB\^-(t>)d(t>. (15) 
By writing dr/dcj) = d{TA + tb) I d(f) — Q , the equation 

rA^(l)A/TA = rB^<i)B/TB (16) 

is easily obtained, similarly to (10), after having written quite naturally = 0^ and 
$ — = ipB- The left side, which contains only terms that arc specific to the first 
segment, is equal to the right side, which contains only terms relative to the second 
segment. We thus exhibit another constant of motion, namely r^d(f)/dT (by shifting 
back to the differential notation), a quantity that turns out to be identified with the 
ratio of the angular momentum L of the particle to its mass m, which we write as 

L/m^r\d(l)/dT). (17) 
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8. Computing the orbit 



Technically speaking in order to determine the trajectory of a moving body free 
from acceleration we apply the following strategy. Knowing the energy E and the 
angular momentum L of the particle of mass m {E and L depend on the initial 
conditions) we can follow the position of that particle by computing the increments of 
its spacetime coordinates t, r and as the proper time r itself advances. Algebraically 
for each increment dr of the proper time we compute (or the computer calculates) the 
corresponding increments dt, dr and d(f) of the coordinate of the mobile body. The 
squares of the increments dt and d4> are extracted from (11) and (17) in the following 
form: 



dt' 



dcP' 



(E/m)2(l-2M/r)-W (18) 
{L/mfr-^dT^ . (19) 

We notice that the expression of dr is missing. We get it by transporting the values 
of dt and d0 into the metric equation (3) and solving it for dr. This yields 

dr' = {{E/rnf - (1 - 2M/r)[l + (L/m)V-']} dr" . (20) 

By dividing both sides of (20) and (19) we directly arrive to the equation of the 
orbit in polar coordinates as 

/ 1 drV _ _ / 2M\ r/m\2 1 



(21) 



9. Trajectory of the planet 



By making the change of variable 

u = 1/r, du = —dr/r' 
equation (21) becomes 



' du 



E' 



m 



= 1 =r,-(i-2M„)l^ + .^ 



(22) 



We now follow the treatment by Weinberg (1972). Consider a test particle moving 
along its closed orbit around the Sun. Its distance r to the central mass M passes 
through a minimum r_ and a maximum r_|_, which correspond respectively to the 
perihelion and the aphelion. In accordance with the change of variable u = (1/r) 
we let 

M_ = l/r_ = f , = l/r_|_ = w , with r_<r<r+,w<u<v. 
At both points the derivative dr/d(f) vanishes, which writes as: 

= (23) 



E' , , m' ' 

j^-{l-2Mv)i — + v' 



w 



0. 



(24) 
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It is easy to extract E'^/ and m^/ from those equations as 

= [v + yj -2M(v^ + VW + w^)]/2M (25) 

— ^ (v + w)(l-2Mv)(l-2Mw)/2M . (26) 
Expression (22) thus takes the form 

^(v + w)(l-2Mv)(l-2Mw)/2M (27) 

\d(j)J 

-(1 - 2Mu)[v + w- 2M(v^ + VW + w^) + 2Mu'^]/2M , 

which does vanish for u = v and u = w. 

Up to now we have not made any approximation but in order to pursue the 
calculation we will take advantage of the smallness of the term Mu and ignore too 
small quantities. In fact the mass of the Sun is 5 x 10~^ s (as was said above) or 
1.5 X 10^ cm and the radius of Mercury's orbit is about 6 x 10^^ cm or 2 x 10^ s, which 
gives a dimensionless factor Mu of order 2.5 x 10~^. Multiplying both sides of (27) by 
the factor (1 — 2Mu)~^ we feel secure to develop that term up to the second order in M 
as 

(1 - 2Mu)~^ ~ 1 + 2Mu + ^M'^u^ . 

We thus get 

(1 - 2Mu)-^ [ ^ j ^{v + w){l- 2Mv){l - 2Mw){l + 2Mu + ^M'^u^)/2M 

-\v + w- 2M{v'^ + VW + w^) + 2Mu'^]/2M . (28) 

The trick to simplify the apparently quite complicated right side consists in noticing 
that we are dealing with a quadratic function of u which vanishes at u = v and u = w 
(if we neglect the terms of order M^). Therefore it has the form 

C{v — u){u — w) . 

The value of the constant C is immediately obtained by letting u — as 
C ^l-2M(v + w), 

which allows us to write equation (28) for the trajectory in the following quite compact 
form 

= (1 - 2Mu) [1 - 2M(v + w)] (v -u)(u-w) . (29) 

\d(j)J 

By taking the square root of both sides of that equation and by neglecting terms 
of order or higher the trajectory of the particle around the Sun can be computed 
by integrating the expression 

# ^ ^ l + M{u + v + w) ^^^^ 

{v — u){u — w) 
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The integration is trivial if one makes the change of variable u ^ ip defined by 

u = -{v + w) -\ — {v — w) COS-;/' (31) 
which easily leads to 

{v — u){u — w) = -{v — w) sin ijj — —du/ dijj 

or 

(/(/ 

^===— = (32) 
J [v — u)[u — w) 

We count the angles starting at the perihelion. At that point = ^ = and u — v. 
Prom that point r increases and thus u — 1/r decreases while the angle -0 increases. 
Taking into account (30), (31) and (32) we arrive at the very simple integral 



m = / 

Jo 



^ 3M. . M. . ■ 
1 H — 7z-[v + w) + —[V — W) COSip 



dip . (33) 



The first term inside the brackets, which is equal to unity, leads to the classical 
Newtonian ellipse. Actually ii (p = ip, expression (31) yields the equation of the 
trajectory in polar coordinates (r, 0) in the form 

u= - — -{v + w) + -{v — w) cos (p , (34) 
7^ 2 2 

which does correspond to an ellipse. Ordinarily the equation is written as 

r = ^ (35) 

1 + e cos (p 

where p is the ellipse parameter (sometimes called the semi-latus rectum) and e its 
eccentricity. 

Identifying (34) and (35) we see that 

11 1 / 1 1 \ 

{v + w) = -{- + -] (36) 
p 2 2 \ r_ r+ / 



and 

V — w r^ — r. 



(37) 



V + w r+ + r_ 
With common notations the major axis of the ellipse is 

2a = r_ + r+ (38) 

and the following relation holds 
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Figure 1. Observation shows that the perihehon of Mercury's orbit advances a little 
with each revolution. This means that the major axis of the ellipse rotates at a tiny 
rate. General relativity accounts for the magnitude of that effect, which is called the 
precession of the perihelion. 



10. Advance of the perihelion 

The second term in integral (33), namely (3M/2)(f + w), reveals that the angle (j) will 
increase more than twice vr when the planet returns to its perihelion, from which it 
started. This phenomenon precisely represents the precession of the orbit. In fact we 
multiply by 2 the increase in the angle between the perihelion r_ (corresponding to 
ijj = 0) and the aphelion r+ (corresponding to ip = it) to obtain a shift per period equal 
to 

GMtt , , 

(39) 



P 

where p is defined in (36). In conventional units we have 
QirG M(g) 



A0 



p (cm) 



The third sinusoidal term in sinip in integral (33) adds only a periodical 
perturbation which produces a kind of noise. 
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Let us calculate the numerical value of the angle of precession. The scmimajor axis 
of the elliptical orbit of Mercure is a = 5.8 x 10^^ cm and its eccentricity is e = 0.206. 
Thus p = 5.55 X 10^^ cm. The mass of the Sun in centimetres is M(cm) = (G/c^)M(g), 
which yields (with M = 2 x lO^''^ g) M = 1.5 x 10^ cm. Thus 

A0 = 5 x 10~^ radian/revolution = 0.103" /revolution 

Knowing that there are 415 revolutions per century we conclude that the advance 
of the perihelion amounts to 

A0 = 43" per century. 
11. Period of revolution 

Let us profit from those formulae as a chance to calculate the period of revolution T 
of a planet around the Sun and recover the Newtonian result. Formula (17) teaches us 
that the increment {l/2)r'^d(f) = {l/2){L/m)dT is proportional to the time increment 
dr. But that quantity is nothing else than the elementary area dA swept out by the 
radius joining the Sun to the planet. Therefore the total area swept out at time r since 
time T = is 

A = (l/2)(L/m)T (40) 

After one complete revolution the total time elapsed is T and the total swept out area 
equals the area of the ellipse, namely Trab if a and b denote the semiaxes. Thus 

7rab= {l/2){L/m)T (41) 

or 

T ^27r(m/L)ab. (42) 
Taking into account (25) which gives (ignoring the M term inside the parentheses) 

= -{v + w)/M^ 



L2 2' pM 
and known relations 

p — a{l — e^) , b = aVl — , 

we find 

r= = ^ . (43) 



M 

In conventional units the formula reads 

"'2 /„3 A„2 



= 47r7(GMg) , 

an expression in which the velocity of light does not appear. This is third Kepler's law. 
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momentum p 




Future path 
of the particle 
(or photon) 



b = impact parameter 



central mass M 




test particle 




with velocity parallel top 

Figure 2. Definition of tlic impact parameter b. The moving particle approaches 
the mass M from a great distance with vector momentum p. A test particle with 
a parallel velocity plunges radially onto the mass M. The distance b between their 
initially parallel paths at "infinity" is the impact parameter b. 



12. Trajectory of light 

The preceding treatment is apparently not relevant to the case of a photon. In fact the 
calculation of the trajectory was done by incrementing the proper time but this latter 
concept has no meaning for a photon since the interval between two events that are 
located on the wordline of a photon is always equal to zero (since at light velocity r = t, 
the interval t"^ =t^ — vanishes). 

Nevertheless it can be shown rigorously (see for instance Taylor and Wheeler, 2000) 
that by letting the mass of the particle tend towards zero, one arrives at the right results. 
Thus for m = our equation (21) takes the form 



That equation will allow us to determine the deviation of light rays passing near the 
Sun. 

It is necessary to specify the parameters found in the formulae. First the angular 
momentum of the moving particle at infinity is equal by definition to the product of 
its linear momentum p by what is called the impact parameter b, which represents the 
distance between the center of attraction (the Sun in the present case) and the initial 
direction of the velocity of the particle (see the figure). 

In other words 




(44) 



L = p b . 



(45) 
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In addition it is known that the momentum p of a photon is equal to its energy E (with 
the units that were chosen). It results at once from this formula that 

L/E = b . (46) 

If the ratio L/E is equal to the impact parameter ft, (44) writes as 

( ]_dry _ ]_ _ f 2M\ 1_ 

[r^dcl)) ~ b'' [ r J ' ^ ^ 

13. Determination of the deflection angle 

Formula (47) will allow us to determine the change in the direction of a light pulse 
caused by the gravitational field of the Sun. To achieve this aim we have to sum up the 
successive infinitesimal increments d0 of the azimuthal angle (p along the path. This 
means that we have to carry out the integration of {d(l)/dr)dr when r varies from the 
minimum distance denoted R {R is the radius of the Sun if the light ray grazes its 
surface) to infinity. We should still multiply that quantity by 2 to account for both 
symmetrical "legs" of the trajectory (the photon first approaches the Sun then recedes 
from it). 

It is necessary to stipulate a further point, namely the relation existing between 
the two quantities b and R that we have introduced, and that are not independent. The 
point r — R corresponds to the place where the light photon is closest to the Sun. There 
the photon moves tangentially. Since at that point there is no radial component, we can 
write that the derivative dr/dt vanishes. It suffices to take the element dr from (47) to 
immediately find 

1 / 2M\ 1 , , 

so that this same equation (47) becomes 

(;^#) " ~ t) '^~{^~ ~) ^ ■ ^^^^ 

The form of the expression dictates to us to pose 
u — R/r 

where u varies between 1 and 0. The last equation (49) then becomes 

{du/d(t)f = (1 - 2M/R) - (1 - 2Mu/R)u^ 
or 

{du/d(l)f = l-u^ - {2M/R){1 - u^) . (50) 

Consequently the infinitesimal variation dcp of the azimuth is given in terms of the 
variation du of R/r hj 

ci0 = [1 - - {2M/R){1 - u^)]-^/^ du 
{l-u^)-y^du 



[1 - {2M/R){1 - m3)(i _ «2)-i]i/2 • 



(51) 
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The presence of the term (1 — u"^) in (51) encourages us to make the change of 
variable 

u — cos a, < u < 1, < a < 7r/2 
which leads to 

d(j)=\l- {2M/R){1 - cos^ a) sin'^ a 

By observing that 

1 — cos^a (1 — cosq;)(1 + cosa + cos^ a) 1 



-1/2 



da . 



(52) 



sin^ a 



— cos a + 



(1 — cos «)(1 + cos a) ' l + cosa 

we end up with the final equation of the trajectory under the form 



1 - {2M/R) (cos a+- 



+ cos a 



-1/2 



da 



(53) 



with 



cosa — R/r . 



It is interesting to emphasize that so far there have been no approximation. This 
is quite rewarding. 



14. Approximations and integration 

The small value of the term M/R will allow us to make an approximation and this 
operation will make us able to complete the integration. In conventional units the mass 
of the Sun is 2 x 10^^ grams and its radius is 7 x 10^° centimetres. By using the factor 
G/(? — l,Ax 10~^^ cm/g which makes it possible to transform grams into centimetres, 
we get the dimensionless quantity 

M/R = (G/c^)M(in grams)/i?(in centimetres) = 2 x 10~^ . 

In (53) we can thus use the classical approximation (1 + e)^ ~ 1 + pe to arrive at 

d(t)= 1 + {M/R) [cos a + Jip^^) ■ (^^) 
Therefore the total variation of the azimuth along the path of the photon is 



I-k/2 



1 + (M/R) (cos a + ^ 

^ ' ^ V 1 + cos a 



da 



M { a 
a + — I sm a + tan — 
R \ 2 



= 2 
= TT + AM/R 



7r/2 



(55) 

(56) 
(57) 



The first term vr gives the total change in the azimuthal angle of the photon where 
there is no Sun present, since in that case the photon follows a straight path. It is the 
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second term that gives the additional angle of deflection A(f) with respect to this straight 
line 

A(j) = AM/R (58) 

or in conventional units 

A(j) = 4(G/c^)M(grams)/i?(centimetres) . (59) 

Numerically at the surface of the Sun (with the values of the mass and the radius 
given above) one flnds A0 = 8, 5 x 10^^ radian, or ( knowing that vr radians equal 180 
degrees and that there are 60 minutes of arc in one degree and 60 seconds of arc in one 
minute of arc) 

A0 = 1.75" . 
15. Summary 

The deflection of light rays by the Sun and the advance of the perihelion of Mercury 
were among the flrst observational conflrmations of General Relativity and are thus of 
major historical value. Unfortunately, in spite of their importance, it is difficult for 
students or non-specialists to access to the full calculations of those effects as they are 
not readily found in the textbooks or in the scientific literature dedicated to Einstein's 
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theory of gravitation. Using a method devised by Taylor and Wheeler (2000) and 
repeating sections of their book, the present work responds to this lack in providing an 
entirely self-contained derivation of the equations and complete algebraic calculations of 
Mercury's precession and starlight deflection. Moreover the paper uses only elementary 
algebra as it does not recourse to tensor formalism. 
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